MPAS Atmosphere
In this recipe, you’ll visualize:
30km MPAS Atmosphere model output
Required Cookbook Sections
TODO
import uxarray as ux
import cartopy.crs as ccrs
grid_path = "../../meshfiles/x1.655362.grid.nc"
data_path = "../../meshfiles/x1.655362.data.nc"
uxds = ux.open_dataset(grid_path, data_path)
uxds
<xarray.UxDataset> Size: 8MB
Dimensions: (Time: 1, n_face: 655362, n_node: 1310720)
Dimensions without coordinates: Time, n_face, n_node
Data variables:
relhum_200hPa (Time, n_face) float32 3MB ...
vorticity_200hPa (Time, n_node) float32 5MB ...<xarray.UxDataset> Size: 8MB
Dimensions: (Time: 1, n_face: 655362, n_node: 1310720)
Dimensions without coordinates: Time, n_face, n_node
Data variables:
relhum_200hPa (Time, n_face) float32 3MB ...
vorticity_200hPa (Time, n_node) float32 5MB ...uxarray.UxDataset
- Time: 1
- n_face: 655362
- n_node: 1310720
- relhum_200hPa(Time, n_face)float32...
- units :
- percent
- long_name :
- Relative humidity vertically interpolated to 200 hPa
[655362 values with dtype=float32]
- vorticity_200hPa(Time, n_node)float32...
- units :
- s^{-1}
- long_name :
- Relative vorticity vertically interpolated to 200 hPa
[1310720 values with dtype=float32]
Show Grid Information
<uxarray.Grid> Original Grid Type: MPAS Grid Dimensions: * n_node: 1310720 * n_face: 655362 * n_max_face_nodes: 10 * n_max_node_faces: 3 * n_nodes_per_face: (655362,) Grid Coordinates (Spherical): * node_lon: (1310720,) * node_lat: (1310720,) * face_lon: (655362,) * face_lat: (655362,) Grid Coordinates (Cartesian): Grid Connectivity Variables: * face_node_connectivity: (655362, 10) * node_face_connectivity: (1310720, 3) Grid Descriptor Variables: * n_nodes_per_face: (655362,)
uxarray.UxDataset.uxgrid
- n_node: 1310720
- n_face: 655362
- n_max_face_nodes: 10
- n_max_node_faces: 3
- node_lon(n_node)float32-175.0 -174.8 ... -103.0 -175.0
- standard_name :
- longitude
- long name :
- Longitude of the corner nodes of each face
- units :
- degrees_east
array([-174.95294, -174.81796, -174.86942, ..., 113.04705, -102.95294, -174.95294], dtype=float32) - node_lat(n_node)float3226.44 26.53 26.67 ... -52.62 -52.62
- standard_name :
- latitude
- long name :
- Latitude of the corner nodes of each face
- units :
- degrees_north
array([ 26.43807 , 26.525747, 26.667755, ..., -10.812317, -52.62263 , -52.62263 ], dtype=float32) - face_lon(n_face)float32-175.0 -103.0 ... -175.2 -174.7
- standard_name :
- longitude
- long name :
- Longitude of the center of each face
- units :
- degrees_east
array([-174.95294, -102.95294, 115.01709, ..., -174.95294, -175.1781 , -174.72784], dtype=float32) - face_lat(n_face)float3226.57 26.57 90.0 ... -52.54 -52.54
- standard_name :
- latitude
- long name :
- Latitude of the center of each face
- units :
- degrees_north
array([ 26.56505 , 26.56505 , 90. , ..., -52.780724, -52.54337 , -52.54337 ], dtype=float32)
- node_face_connectivity(n_node, n_max_node_faces)int640 163843 163844 ... 655360 655361
- cf_role :
- node_face_connectivity
- long name :
- Faces that neighbor each node
- start_index :
- 0
- _FillValue :
- -9223372036854775808
array([[ 0, 163843, 163844], [ 0, 163844, 163845], [ 0, 163845, 163846], ..., [655353, 655354, 655355], [655356, 655357, 655358], [655359, 655360, 655361]]) - face_node_connectivity(n_face, n_max_face_nodes)int640 1 ... -9223372036854775808
- cf_role :
- face_node_connectivity
- long name :
- Maps every face to its corner nodes.
- start_index :
- 0
- _FillValue :
- -9223372036854775808
array([[ 0, 1, 2, ..., -9223372036854775808, -9223372036854775808, -9223372036854775808], [ 5, 6, 7, ..., -9223372036854775808, -9223372036854775808, -9223372036854775808], [ 10, 11, 12, ..., -9223372036854775808, -9223372036854775808, -9223372036854775808], ..., [ 983022, 1228798, 983034, ..., -9223372036854775808, -9223372036854775808, -9223372036854775808], [ 1310719, 983028, 983033, ..., -9223372036854775808, -9223372036854775808, -9223372036854775808], [ 1310719, 983039, 983038, ..., -9223372036854775808, -9223372036854775808, -9223372036854775808]])
- model_name :
- mpas
- core_name :
- init_atmosphere
- source :
- MPAS
- Conventions :
- MPAS
- git_version :
- v6.0-dirty
- on_a_sphere :
- YES
- sphere_radius :
- 6371229.0
- is_periodic :
- NO
- x_period :
- 0.0
- y_period :
- 0.0
- history :
- Mon Aug 16 12:09:52 2021: ncks -v latVertex,lonVertex,latCell,lonCell,latEdge,lonEdge,cellsOnVertex,verticesOnCell,verticesOnEdge,edgesOnVertex,edgesOnCell,nEdgesOnCell,angleEdge x1.655362.init.nc foo.nc mpirun -n 18 ./init_atmosphere_model
- parent_id :
- 0hhmjs26t0
- mesh_spec :
- 0.0
- config_init_case :
- 7
- config_calendar_type :
- gregorian
- config_start_time :
- 2016-08-01_00:00:00
- config_stop_time :
- 2016-09-10_00:00:00
- config_theta_adv_order :
- 3
- config_coef_3rd_order :
- 0.25
- config_num_halos :
- 2
- config_nvertlevels :
- 75
- config_nsoillevels :
- 4
- config_nfglevels :
- 38
- config_nfgsoillevels :
- 4
- config_months :
- 12
- config_geog_data_path :
- /glade/p/work/wrfhelp/WPS_GEOG/
- config_met_prefix :
- FILE
- config_sfc_prefix :
- SST
- config_fg_interval :
- 86400
- config_landuse_data :
- USGS
- config_topo_data :
- GMTED2010
- config_use_spechumd :
- NO
- config_ztop :
- 40000.0
- config_nsmterrain :
- 1
- config_smooth_surfaces :
- YES
- config_dzmin :
- 0.3
- config_nsm :
- 30
- config_tc_vertical_grid :
- YES
- config_extrap_airtemp :
- linear
- config_static_interp :
- NO
- config_native_gwd_static :
- NO
- config_gwd_cell_scaling :
- 1.0
- config_vertical_grid :
- YES
- config_met_interp :
- YES
- config_input_sst :
- NO
- config_frac_seaice :
- YES
- config_pio_num_iotasks :
- 0
- config_pio_stride :
- 1
- config_block_decomp_file_prefix :
- x1.655362.graph.info.part.
- config_number_of_blocks :
- 0
- config_explicit_proc_decomp :
- NO
- config_proc_decomp_file_prefix :
- graph.info.part.
- file_id :
- ytj0fh0qmd
- NCO :
- netCDF Operators version 4.9.5 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)
Relative Humidity
uxds['relhum_200hPa'][0].plot(projection=ccrs.Robinson(), pixel_ratio=4.0, features=['coastline'])